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Wc consider the transport properties in an applied magnetic field of the spin 5 = 1/2 Kondo 
model with ferromagnetic exchange coupling to electronic reservoirs, a description relevant for the 
strong coupling limit of underscrcened spin 5 = 1 Kondo impurities. Because the ferromagnetic 
Kondo interaction is marginally irrelevant, perturbative methods should prove accurate down to 
low energies. For the purpose of this study, we use a combination of Majorana diagrammatic 
theory with Density Matrix Numerical Renormalization Group simulations. In the standard case of 
antiferromagnetic Kondo exchange, we first show that our technique recovers previously obtained 
results for the 7"-matrix and spin relaxation at weak coupling (above the Kondo temperature). 
Considering then the ferromagnetic case, we demonstrate how the low-energy Kondo anomaly splits 
for arbitrary small values of the Zeeman energy, in contrast to fully screened Kondo impurities near 
the strong coupling Fermi liquid fixed point, and in agreement with recent experimental findings for 
spin S = 1 molecular quantum dots. 



I. INTRODUCTION 

Magnetotransport measurements in quantum dots al- 
low to probe and tune the magnetic properties of nanos- 
tructurcs in a situation where the conductance is the only 
experimentally accessible quantity Cotunneling Zee- 
man spectroscopy under an applied voltage bias {i.e. in 
a process that docs not electrically charge the quantum 
dot) can be used for instance to resolve the magnetic 
degeneracies of quantized energy levels, in similar spirit 
to usual optical analysis of atomic and molecular spectra 
[see Refs.^'^ for an illustration in a two-electron quantum 
dot showing nearly degenerate singlet and triplet states] . 
A striking difference between quantum electronic trans- 
port and quantum optics is the easiness to realize the 
strong coupling of a discrete quantum system to a macro- 
scopic number of external degrees of freedom, simply by 
contacting a single atom or molecule to metallic elec- 
trodes. While standard physical effects such as screening 
of the molecular charging energy^ and rcnormalizations 
of the atomic energy levels from quantum tunneling^'^'^ 
can occur, these electronic setups opens a window to re- 
alize even more subtle many-body physics. One well- 
studied example is the Kondo effect in strongly tunnel- 
coupled quantum dots, where the magnetic screening of a 
spin active discrete level is performed by the bath of sur- 
rounding conduction electrons*^^^. Interestingly, a va- 
riety of different and possibly exotic Kondo effects can 
occur in simple extensions of the Kondo model, several 
of which have been already realized experimentally [see 
Ref. 4 for a review]. While the standard Kondo effect 
involves a single spin 5 = 1/2 coupled to a single elec- 
tronic bath and leads to Fermi-liquid properties at low 
temperatures, a variant as simple as the spin S = 1 case 
already presents anomalous transport properties"'^^'''^^. In 
the latter situation, conduction electrons cannot manage 
a full compensation of the magnetic impurity, resulting 
in a partially screened ground state with remanent log(2) 
entropy^^"^'', yet unitary conductance at zero tempera- 



ture G = 2e^//i (e is the electron charge and h Planck's 
constant) . At low energy, the remanent magnetic degrees 
of freedom decouple from the electronic reservoirs (in a 
logarithmically slow fashion) and can be recovered by the 
application of an arbitrarily small Zeeman energy In 
contrast, the fully screened case of a spin S = 1/2 leads 
to a singlet ground state with small but non-zero bind- 
ing energy (the so-called Kondo temperature), which re- 
mains stable to the application of a magnetic fields that 
do not exceed the Kondo energy. These simple exam- 
ples illustrate nicely how Zeeman spectroscopy can be 
used to distinguish between different types of many-body 
states^^'^*. 

The purpose of this paper is to investigate in detail 
the magnetotransport properties of such highly spin- 
polarizable underscrcened 5 = 1 Kondo impurities. We 
will make use here of the Nozieres-Blandin description 
of the strong coupling fixed point in terms of the spin 
5=1/2 Kondo model with ferromagnetic Kondo ex- 
change^^'^^. In that case, the magnetic interaction be- 
tween the quantum impurity and the reservoirs can be 
treated by perturbative methods, allowing a controlled 
description of the transport properties. While one could 
use here the functional renormalization group method- 
ologies previously developed to deal with the weak cou- 
pling regime of the antiferromagnetic case^^~^^, we will 
show that a comparable description can be achieved by 
more standard diagrammatic methods, thanks to an use- 
ful mapping of a single spin S = 1/2 in terms of Majorana 
fermions^'^"^^. The advantage of the Majorana technique 
with respect to more usual representations of the spin al- 
gebra lies in the recent understanding that spin dynamics 
can be described by one-particle Majorana Green's func- 
tions, allowing the use of Dyson's equation within the 
standard Feynman diagram expansion^^'^^'^^. We will 
demonstrate that transport properties, as described by 
the T-matrix of the Kondo problem, become similarly 
much simpler in the Majorana language. Indeed, while 
standard spin representations would require the compu- 
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tation of six-point correlation functions, the Majorana 
description expresses the T-matrix as a four-point (two- 
body) fermionic correlator, where the RPA expansion for 
response functions can be used at leading logarithmic or- 
der. In order to test the method, we will first discuss how 
the usual weak-coupling description of the antiferromag- 
netic Kondo problem can be done within the Majorana 
diagrammatics, both at zero and finite magnetic field. 
A comparison to calculations based on the full density 
matrix extension^^''^'' of the Numerical Renormalization 
Group (NRG)^^"^^, and previous analytical results^^'^^ 
will be made. We will then compute in detail the spin 
and transport dynamics of a ferromagnetically coupled 
quantum dot, both within the Majorana description and 
the NRG simulations, which are shown to match quan- 
titatively, apart from a narrow domain near the Zeeman 
energy, where the NRG data lack sufficient energy reso- 
lution. Our main result concerns the Zeeman splitting of 
the T-matrix in the ferromagnetic case, which we show 
to occur for arbitrarily small values of the applied mag- 
netic field, in agreement with recent experimental find- 
ings for underscreened S = 1 Kondo impurities in molec- 
ular quantum dots'^'^. We also note that the ferromag- 
netic Kondo model considered here has also been recently 
proposed for some types of atomic size metallic constric- 
tions on the basis of ab-initio calculations'^'*. 

The paper is organized as follows. In Sec. II we in- 
troduce the Kondo model relevant for transport in quan- 
tum dots, discussing transport and magnetic properties. 
The differences between fully and partially screened mo- 
ments, with the connection to the ferromagnetic Kondo 
problem, will be presented. In Sec. Ill, we develop the 
general Majorana diagrammatic method, both at zero 
and finite magnetic field. Results arc finally given in 
Sec. IV, where both antiferromagnetically and ferromag- 
netically coupled spin 5=1/2 impurities are considered, 
using the Majorana technique and NRG calculations. A 
critical discussion of numerical and analytical scheme is 
given in the conclusion. 



II. KONDO MODELS FOR TRANSPORT IN 
QUANTUM DOTS 

A. Spin S = 1/2 case 

The conductance Q oi a quantum dot is the main phys- 
ical property that can be directly measured in quantum 
dot experiments. In semiconducting dots, electrons are 
confined to a small region with the help of top-gates, lead- 
ing to a set of discrete energy levels subject to electron- 
electron interaction effects, and tunnel coupled to macro- 
scopic reservoirs. Under a voltage bias, a finite current 
flows through the dot, and the condiictance is deter- 
mined by the transition rate of tunneling events from 
source to drain. The concept of single-electron transis- 
tor comes from the capacitive coupling in small confined 
dots, where the discrete number of electrons on the dot 



can be tuned using a backgate-voltagc-'^'^. If the number 
of trapped electrons on the dot is odd, an electronic con- 
figuration with an unpaired spin obtains, leading usually 
to a S* = 1/2 state. The tunneling into the reservoirs then 
results in an effective antiferromagnetic exchange inter- 
action J > 0, which tends to suppress the magnetic state 
of the dot. The physical origin of the antiferromagnetic 
exchange is Pauli's principle, which favors second order 
tunnel processes for antiparallel spin orientation of the 
electron in the dot and the ones in the reservoirs. This 
interaction can be described by the Kondo Hamiltonian: 

Qia2 ki k2 <^i<^2 

where r denotes the vector r — (t^, t^, t^) formed by 
the three Pauli matrices. For simplicity we have assumed 
here symmetric coupling to the left and right reservoirs. 
Besides the quantized spin S in the dot, we have in- 
troduced fermion operators describing electrons in the 
reservoirs, as given by the Hamiltonian: 

^0 = X]^feCafeaCafea' (2) 
kaa 

where the operator cj^,-.^ creates one electron in the lead 

a = L, R with momentum k, energy and spin a =t, i- 
It is useful to define a "local" electronic state of the 
reservoirs cj^ = J2koc4ika/'^^ allowing to write the 
Kondo interaction in a more compact way: 

Hk = JS .^cL.^Co,,. (3) 

In this form, a single screening channel couples to the 
impurity spin, which allows to understand the ultimate 
formation of a non-degenerate singlet ground state. The 
conductance through the dot can be evaluated in the lin- 
ear response regime using Kubo formula^^'^^ 

1 r°° 

a=lim-/ e-*([/(t),J(0)])cit, (4) 
with the current operator: 

in terms of the charge operator in each lead. The 
complete scattering process can be formally described 
by the full Green's function of the conduction elec- 
trons" 

G^^,{k,k\u) = Gg{k,oj)dk,k' (6) 
+Gg{k,co)rJ'{co)Gl{k\u>), 

where TJ^{lo) denotes the retarded T-matrix (which is 
momentum independent due to the local nature of the 
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Kondo interaction). We have also introduced the free 
Green's function GQ^{k,oj) of the conduction electrons: 



(7) 



In all generality, the equilibrium conductance can be re- 
lated to the T-matrix by the formula (here in the case of 
left /right symmetric barriers): 



2e2 
h 



with the density of states at the Fermi level po = 

— (l/7r)XT77. 1 / (— gfe+zO'^), and the Fermi function np- 
In what follow, the standard notation T ^ = Im[T^] is 
used. Now focusing on the Kondo problem, straightfor- 
ward derivation using equation of motion for the conduc- 
tion electron Green's function (6) provides the explicit 
form of the T-matrix: 



J 



{S^} 



E 



(9) 



X (Tr 



5(0) 



using standard imaginary time-ordered correlations func- 
tions. The T-matrix in Eq. (9) describes spin-flip pro- 
cesses to all orders, but its pcrturbative expansion in the 
small parameter pqJ is known to be logarithmically di- 
vergent, giving rise to the so-called Kondo problem. The 
mathematical basis for the systematic resummation of 
logarithmic terms in perturbation theory lies in renor- 
malization group ideas, where electron states of the reser- 
voirs are skimmed progressively from high to low energy. 
At leading order, this procedure is encapsulated by the 
following flow equation^'^^: 



(10) 



where D is the running high frequency cutoff of the con- 
duction band, starting with the initial value I^o at the 
band edge (so that po = l/(2-Do) for a purely flat density 
of states). Direct integration leads to the renormalized 
dimensionless Kondo interaction at the energy scale D: 



PoJr 



PqJ 
1 - poJlog 



1 



Do 
D 



log 



D 
Tk 



(11) 



with the Kondo temperature Tk = Doe~poJ which sets 
the frontier beyond which perturbation theory breaks 
down. In the weak-coupling regime T ^ Tk the con- 
ductance Q{T) is given by: 



g{T) = 



2e2 Stt^ 



1 



h 



16 log^ 



(12) 



The divergence at temperature of the order of Tk marks 
the formation of a Kondo resonance, that in reality satu- 
rates in the low temperature strong- coupling regime to- 
wards the unitary limit Capturing the crossover 

from high to low temperature is possible only by non- 
perturbative methods, such as Wilson's numerical renor- 
malization group (NRG)'^^. The resulting physical pic- 
ture of the strong-coupling regime is the screening of the 
local spin S = 1/2 by the surrounding conduction elec- 
trons^^, with a binding energy for the singlet ground state 
of the order of Tk- Scattering in the low-temperature 
regime is then dominated by irrelevant terms, leading to 
a Fermi- liquid behavior for T <^ Tk- 



g{T) 



T 

16 \Tk 



(13) 



Alternatively to the temperature dependence of the equi- 
librium conductance, magnetic field and finite bias effect 
lead similarly to a suppression of transport when the scale 
of the perturbation exceeds Tk- Focusing on the finite 
bias Kondo resonance (at temperatures much below Tr-), 
it has been established theoretically and experimentally 
that the Kondo anomaly is robust to a Zeeman field of 
the order of Tk, which reflects the singlet binding of the 
groimd state. One goal of the paper is to study the Zee- 
man splitting of Kondo resonances associated to a spin- 
degenerate ground state. A basic model for the formation 
of this more exotic state is reviewed now. 



B. Spin S = \ CEise 

The Kondo Hamiltonian (1) is straightforwardly gener- 
alized to higher spin values. In particular, quantum dots 
with an even number of electrons may lead to the forma- 
tion of a spin S = \ magnetic state, which can undergo 
the so-called underscreeed Kondo effect^'^'^'^^"^''. In the 
weak coupling regime [See Ref. 22 for a thorough recent 
study], renormalization of the antiferromagnetic Kondo 
interaction still follows Eq. (10), and the conductance 
reads for T > Trt: 



Q{T) 



2e2 1 



h 8 



log' 



(14) 



Thus a Kondo resonance develops at T < Tk which can 
reach up to 2e^/h values at low temperature^'^. The 
crucial difference with respect to the case S = 1/2 is 
that a spin S = 1 coupled to a single screening chan- 
nel cannot bind into a total singlet, so that the result- 
ing ground state has to accomodate a log(2) entropy. 
Nozieres and Blandin^^ showed early on that the scatter- 
ing process on such a partially screened state was anoma- 
lous. Their argumentation was based to a description 
of the strong coupling fixed point in terms of an effec- 
tive spin S'^^ = 1/2 describing the remaining degeneracy 
in the system. Because the original spin S = 1 tends 
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to bind antiparallel spin electron states, Pauli principle 
leads to fluctuations with parallel spin electrons, thus 
leading to an effective ferromagnetic interaction < 0, 
with (X —[PoJ]~^ in the mathematical limit poJ 1 
considered by Nozieres and Blandin^^. Integrating the 
flow equation (10) with a ferromagnetic coupling leads 
to the renormalized Kendo interaction: 



PoJ' 



eff 



l + Po|J^«^|log 



Do 
D 



(15) 



which does not diverge upon decreasing the cutoff D, but 
rather logarithmically vanishes at low energy. Consider- 
ing the lowest order contribution to transport, one then 
arrives at the following result for the conductance in the 
strong coupling regime T <C Tk- 



g{T) 



2e^ 



1 - 



37r^ 
16 



1 



log' 



T 
Tk 



(16) 



This slow logarithmic approach to the low temperature 
fixed point was observed in recent experiments^^' Our 
goal for the rest of the paper is to understand the magne- 
totransport of such underscreened Kondo anomalies. In- 
deed, despite the binding energy Tk, the underscreened 
ground state is still two-fold ciegenerate, and is expected 
to be highly sensitive to the Zeeman effect (the spin sus- 
ceptibility is strictly speaking infinite at zero tempera- 
ture). Before addressing this physical question, we will 
develop the required methodology. 



III. MAJORANA DIAGRAMMATICS FOR THE 
PERTURBATIVE KONDO PROBLEM 

A. Majorana representation and observables 

Calculating the perturbative expansion in the Kondo 
problem (1), especially at high orders, requires the eval- 
uation of correlation functions with multiple occurence 
of the spin operators. These involve cumbersome algebra 
due to the peculiar commutation relations of spin op- 
erators (which are neither fermions nor bosons), which 
naively precludes the use of the useful Wick's theorem. 
A standard alternative is to represent the quantized spin 
with fictitious particles, for instance using the Abrikosov 
fermion^'^ description of a spin S = 1/2: 



(17) 



with a spin-dependent single fermion level submitted 

to the constraint flfa = 1- While this now allows the 
use of Wick's theorem, the computation of observables is 
not very much simplified, because a single spin opera- 
tor involves the product of two fermions. For instance 
the spin correlation functions Xii^) = (with 



i = x,y,z) involve four fermions, while the 7~-matrix 
Eq. (9) requires a correlator with six fermions (includ- 
ing the two electrons originating from the leads). Order 
by order terms can thus be obtained with increasing ef- 
fort, but systematic ways to perform resummations of the 
perturbation series seems hopeless, especially for the T- 
matrix, so that the perturbative renormalization group is 
often preferred in practice. Despite the successes of the 
renormalization ideas, difficulties arise when considering 
physical quantities which depend on several energy scales 
together (such as temperature, magnetic field, voltage, or 
frequency), in which case a single-scale renormalization 
procedure does not apply, and more elaborate functional 
forms of the RG must be used^'^"^^. While these method 
will certainly perform well for the ferromagnetic Kondo 
problem that we consider here, we would like to consider 
an alternative methodology, which, while not as general 
and powerful as the RG, will lead to similarly accurate 
results by the much simpler means of standard Green's 
function perturbation theory. 

The idea that we will follow, building on earlier 
works^^"^'' is to use a Majorana fermion representation 
of a single spin S = 1/2: 



T»7 X "n, 



(18) 



where r} = (??i, »72, %) is a triplet of real fermions which 
satisfy the anticommutation relation {T?a,%} = 5ab- One 
first advantage is for instance that the Hilbert space does 
not need to be restricted, because 5' = | is automati- 
cally fulfilled^^"^^. Despite the writing of the spin oper- 
ators in terms of a product of two Majorana fermions, a 
real gain can be achieved by the alternative decomposi- 

tion24.25. 



S = $T7 with $ = — 2i?7i?72?73. 



(19) 



One can then easily check that the fermion $ commutes 
with the Kondo Hamiltonian, so that for instance the 
^-axis spin correlation function can be written: 



X.(r,0) = -(r,[5^(r)5^(0)]) 



$(T)773(r)$t(o)r,t(o) 
= G$(T)G^3(r) 
introducing the Majorana Green's functions: 



(20) 



(21) 



as well as the propagator of the fermion which remains 
free to all orders in perturbation theory: 



G^{t) = -{Tr [$(r)$t(0)]) = -^Sign(T). 



(22) 



Now the spin susceptibility Eq. (20) is captured just by 

a single Majorana propagator, which allows the use of 
Dyson's equation instead of more cumbersome response 
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FIG. 1. Interaction vertex associated with the first term in 
the Kondo Hamiltonian (24) expressed in the Majorana lan- 
guage. Full lines denote conduction electrons and broken lines 
Majorana fermions. 



function framework. Although Majorana fermions are 
real wc distinguish formally r]^ and rj. This step makes 
the diagrammatic expansion clearer once expressed in 
terms of particle and holes^^. 

Considering now the T-matrix of the Kondo model, 
originally a six-fermion correlator, the above arguments 
leads to a simpler four-fermion function: 



(23) 



<7iO'2 ah 



X (T.. 



Up to the constant $ propagator, the T'-matrix assumes 
the form of a response function, for which standard per- 
turbative methodology exists, although more involved 
than for the Majorana propagators required for the com- 
putation of spin dynamics. 

In order to proceed with perturbation theory, we finally 
need the Majorana form of the Kondo Hamiltonian (1): 



H 



K 



-iJ 



(24) 



-iJ 



Interaction vertices between pairs of Majoranas and pairs 
of conduction electrons are readily identified from this 
rewriting of the Hamiltonian, as shown in Fig. 1. 



B. Majorana diagrammatics at zero magnetic field 

In this subsection we will consider the perturbative ex- 
pansion of the T-matrix in terms of Majorana fermions 
given by Eq. (23). The required response functions ap- 
pearing in the T-matrix are given diagrammatically in 
Fig. 2 where the box represents spin-spin interactions 
to all orders. Since there are three different Majorana 
fermions, in principle nine different correlators have to 
be expanded in perturbation theory, but in reality only a 
small number of independent correlation functions arise, 
as we will see. To show this, we define the following ab- 




FIG. 2. Graphical representation of the general correlator 

involved in the Majorana decomposition of the T-matrix 
Eq. (23) and Eq. (25). The shaded box indicates interac- 
tion processes between conduction electrons and Majorana 
fermions to all orders in perturbation theory. 



breviation for the response functions introduced in Fig. 2: 

/ iivii vi) *(i vi i m) (i m t '?3)\ 
m(t) = r]2 i m) a m i m) -«U m t , 
V {tmim) i{\v3iv2) (t??3t??3)/ 

which allow to write more compactly the T-matrix: 



Ttir) 



'V 



G$(r)(l, 1, 1) •M(r)- 1 



(26) 



,1, 



In order to perform a leading-log calculation of the T- 

matrix, we will use in what follows a RPA resummation 
with a single polarization bubble, owing to the fact that 
non-RPA diagrams (i.e. involving crossings) are higher 
orders in logarithmic singularities. In Fig. 3 this RPA ex- 
pansion is shown for the response function entering the 
T-matrix. In the case of spin-rotation invariance, all the 
bubbles in Fig. 3 take the same value (given here in imag- 
inary time): 

lionir) = Goc(T)Go^(r). (27) 
Here the bare Majorana Green's function simply read: 



Got,(«w) 



1 



(28) 



while the conduction electron propagator are given by: 

1 



Goc(iw) = E 



(29) 



The first terms in the diagrammatic series for the T- 
matrix are easily checked to read: 



Tt(T) = ^G^(r) 



^Hor, + ^6Ho,, * Uor, + ... 



(30) 

The factors 3 and 6 stem from equation (26), because at 
order the matrix M in Eq. (25) is purely diagonal. 
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FIG. 3. RPA-expansion of the response functions entering the T-matrix in the Majorana language. The single polarization 
bubble involves the product of a free Majorana fermion propagator Goi; and of a conduction electron propagator Goc, and the 
RPA resummation corresponds to the leading logarithmic approximation of the standard one-loop renormalization procedure. 



while at order it is purely off-diagonal. The convolu- 
tion operator in Eq. (30) is standardly defined as: 



dri Uo{t - Ti)no,, 
Jo 



(ri) (31) 



We can already anticipate that a RPA-like resummation 
will be possible by Fourier transforming to imaginary fre- 
quencies. However, the numerical coefficient of the se- 
ries (30) need to be obtained, and for this we will use a 
recursion based on the following diagrammatic rules. In 
principle, there are only two different response-functions. 
The first type is a response-function where the same Ma- 
jorana fermion enters and leaves the interacting region 
(gray box of the T-matrix in Fig. 3), or equivalently when 
no spin-flip of the conduction electrons occurs. These 
contributions are the diagonal terms of the matrix given 
in equation (25) and we call a„ their combinatorial factor, 
so that they contribute in total to 3a„ to the T-matrix in 
Eq. (26). The second type is a response- function where 
different Majorana fermions enter and leave the interact- 
ing region, so that the conduction electrons have flipped 
their spin. These processes are labeled with a combi- 
natorial factor bn associated to the non-diagonal terms 
of the matrix in Eq. (25), thus contributing 66„ to the 
T-matrix, so that we have the general RPA-series: 



q j2 

rt=0 



26„) (^-j [Ho,, Ho,,] 

(n-j-l) times 

(32) 

where oq = 1, 6o = 0, ai = and hi — 1. In order 
to find the explicit expression for a„ et we note the 
combinatorics of Fig. 4 which follows from the fact that 
a Majorana fermion cannot scatter into itself, due to the 
identity 77^ = 1/2. Thus a diagonal process at order n 
(described by a„) results necessarily from a combination 
of all non-diagonal processes at order n — 1 (described 
by bn-i) and of a single off-diagonal term, for which two 
possibilities arise. This shows readily that a„ 25„_i. 
Similarly an off-diagonal process at order n can be de- 
composed either by an arbitrary diagonal processes at 
order n — \ and a single off-diagonal one, or by all possible 
non-diagonal processes at order 71 — 1 followed by a sin- 
gle non-diagonal one, so that 6„ = a„_i + 6„_i. We can 
trivially solve this recursion, and we find a„ + 26„ — 2". 
Using this result, the T-matrix is given at RPA level by 



K 



0® 
®® 

dXD 
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a„ — 26„_i 



bn — a-n-l + bn-l 



FIG. 4. (color online) Combinatorics of the RPA diagrams 
at zero magnetic field: circles denote a single polarization 
bubble involving a Majorana fermion with flavor 1,2, or 3. 
The coefficient a„ counts the number of possibilities at order n 
to end up with the same index, while the coefficient 6„ counts 
the sequences with different head and tail. The mathematical 
recursion on the right side follows from this obvious pictorial 
description. 



q t2 

Tt(r) = _G*(r)n,(r), 



(33) 



where n^(r) is simply determined by a geometric series 
in Matsubara frequencies: 



no,,(w) 

ILjUJ^ = ; r. 

' i-jno^(zi^) 



(34) 



We finally compute the expression for the T-matrix, 
first by Fourier transforming to fermionic Matsubara fre- 
quencies: 



3J2 ^ 

T\{i(jJn) = — -^G'$(wi)n(ia;„- wi) 



3J_ 
Stt 



2 poo 



dw2 n/(w2) 



l + 2nB(-W2) 



(35) 



W2 - ^^^n 



and using a spectral decomposition (ng is the Bose dis- 
tribution). The analytic continuation is then straight- 
forward, and we obtain the imaginary part of the real- 
frequency T-matrix: 



3j2 

T ^(w) = coth 



/3w 



Im 



(36) 



with (3 = 1/T the inverse temperature. The polarization 
diagram Ho,, is computed by similar means from Eq. (27) 
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FIG. 5. (color online) T-matrix of the Kondo model with 
antiferromagnetic exchange at zero magnetic field, comparing 
the NRG simulations (full line) to the poor man's scaling one- 
loop result (dotted line) and to the Majorana diagrammatics 
at RPA level (dashed line) given by Eq. (38). 



and reads for a flat conduction band with density of states 
Po ~ l/(2£') and half-bandwidth D: 



n«(c.) = ^iog 



-fi7r^Sign(w)e(w2_^2) 
(37) 



leading to the final result for T-matrix at T = 0: 



3^ 
16 



1 - Y log 



(po J)^e(c.^ - D^) 
1 



I o^- 



4 



(38) 

Inserting this expression into the conductance for- 
mula (8), we recover the standard high-temperature re- 
sult Eq. (12) for T ^ Tk- Let us contrast our ex- 
pression to the standard weak-coupling renormalization 
group result in the case of antiferromagnetic Kondo ex- 
change: —t:pqT-^'^{lu) = '^[log{uj/Tx)]~'^ ■ In the range 
Tk ^ a; <C -D, both are completely equivalent. For fre- 
quencies near the high energy cutoff D, the correct band- 
edge contribution appears in Eq. (38), which is usually 
neglected in the RG flow (but can be in principle in- 
corporated). Stronger differences occur in the regime 
near Tk, where a maximum occurs in Eq. (38) instead 
of a divergence. The reason is the presence of a con- 
stant offset of order in the denominator, which cuts 
off the Kondo divergence. This term would result from 
next-to-leading order contributions at two-loop order in 
the standard RG. We emphasize that perturbation still 
breaks down in Eq. (38) because the effective dimension- 
less Kondo coupling has reached values of order 1. Yet, 
the convergence of our result to NRG data seems slightly 
better than the lowest order RG expression, see Fig. 5. 
We now show how the Majorana diagrammatics can be 
extended to include the Zeeman effect as well. 



C. Majorana diagrammatics at finite magnetic field 

The Zeeman effect can be easily studied in quantum 
dots, by applying a magnetic field in the plane in the case 
of semiconducting two-dimensional electron gas devices. 
We note that the latter systems have relatively small g- 
factor, but because of the low Kondo temperature that 
they can achieve (in the range lOOmK to IK), a compe- 
tition with the Kondo effect can occur at magnetic fields 
of a few Tesla. In molecular quantum dots larger g ~ 2 
values are obtained, which helps realizing the competi- 
tion between Kondo and Zeeman effects. In this section, 
we simply add the Zeeman term to the Kondo Hamilto- 
nian (1) 



(39) 



where B is measured in units of the Zeeman splitting. We 
note that the Majorana fermion propagator Gab{T) has 
now non-diagonal elements, which makes it more difficult 
to carry out the previous RPA for the T-matrix. The 
first idea is to transform the Majorana fermions so that 
their propagator becomes diagonal, which can be done 
by introducing one real fermion: 

The Kondo Hamiltonian is given after this transforma- 
tion by: 

'^(^/Vco^-^/Vcoil (41) 



{f\f] = i- 



(40) 



2 

The last term, corresponding to a small magnetic field 
acting on the bath of conduction electrons, is expected 
to remain perturbatively small. However, we will see in 
Sec. IV C that it needs to be kept in order to describe 
properly the self-energy in second order. The Zeeman 
term now becomes purely diagonal: 

Bff - \b, (42) 



where the last term is only a constant and will be dropped 
from now on. Thus the transformation leads to two dif- 
ferent fermionic propagators 

Gf[iuJn)^-. — G,,3(icj„) = — , (43) 

where G/ is a canonical fermion propagator and Gn^^ is 
the propagator of the unchanged Majorana fermion 773. 
The label 3 can be dropped and we will use in what 
follows G,, as the Majorana fermion propagator. Next, 
the T-matrix is expressed in terms / and 77: 

rt(T) = - ^g$(t) (2a / ; /t) + \/2a / 1 v) 



+ V2(tr,;/t) + (tr7try)), 



(44) 
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FIG. 6. (color online) Lowest order contributions to the RPA 
evaluation of the T-matrix in the case of finite magnetic field. 



We have now two different bare polarization bubbles Ho / 
and Ho,,, associated to / and rj respectively, which are 
readily evaluated: 



noj/t^) = y log 
no/(w) = Polog 



D + B -uj 



+ i^Trpo Sign(iLj) (45) 



nrpo0{LJ-B), (46) 



in the case B > 0. Because of these two different con- 
tributions, the combinatorics cannot be guessed as easily 
as in the previous section, so we must use a more general 
resummation technique. The general structure of per- 
turbation theory (with bare polarization bubbles) can be 
considered at the light of Fig. 6. We introduce a dimen- 
sionful coefficient An which encodes processes at order 
n that start and finish with the bubble Ho/, and a co- 
efficient Bn that involve Hoy and Ho,, at its extremities. 
Note that these coefficients are not only combinatoric in 
character (like the coefficients a„ and 6„ defined in the 
previous section), but also include a sum of terms of the 
type [no/]P[no,,]"~'' with integer p, due to the fact that 
we have two independent polarization bubbles now. The 
correlations appearing in Eq. (44) are then expressed as: 



itvtv)^ Ho,, + 2 



'J 



2 oo / 7\ 

n=0 ^ ^ 



V 2 

OO ^ 

V2 ((t 7] iP) + {ift V)) = 2V2j2Bn 



(47) 
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A„ =no/(%/2S„_i + A„_i) 



S«=®® ® Bn = V2UoAn^, 

FIG. 7. (color online) Combinatorics of the RPA diagrams at 
finite magnetic field: circles denote a single polarization bub- 
ble involving a real fermion / or a Majorana ry. The coefficient 
An picks the contributions at order n that end up with the 
same index /, while the coefficient collects the sequences 
with different head and tail. The mathematical recursion on 
the right side follows from this obvious pictorial description. 



SO that the complete RPA series (44) for the T-matrix 
reads: 

T;'^H = -yCoth(^)n"«(c.) 



n^iuj) =nS,(..) + 2£(0"|A: 

+ V2Bn I. 



An 



(48) 



A recursion is easily established on the basis of Fig. 7, 
which can be best cast and solved in matrix form: 



Ho/ 


\/2no/\ 


[A . 


V2U0r^ 


) 




no/ 













(49) 



with Aq = Ho/ and -Bg = as can be inferred from 
the lowest order perturbative expansion shown in Fig. 6, 
Inserting Eq. (49) into the expression (48) for the T- 
matrix, we end up with the final compact result at zero 
temperature: 



t;«(^) 



J2 



Sign(a;) Im 



n, 



077 



SHq/ -)- I JIIo/IIo^ 



(50) 

which will be investigated in Sec. IV. If we take the 
limit i? — >■ before the limit of zero temperature, we find 
Ho,, = Ho/, so that expression (50) reduces to 



T;«(w) = -ySign(c.)Im 
= — — Sign(a;) Im 



3no^ (1 + jUp,,) 

(1 - jHo^) (1 + |no„) 

SIIo,, 



(51) 



which is exactly the one established previously in Eq. (36) 
for zero magnetic field. Note however that the limit 
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B and T — >■ clearly do not commute by comparing 
Eq. (45-46) and Eq. (37), due to the abrupt change in 
the ground state degeneracy of the free impurity spin as 
soon as a small magnetic field is turned on. At zero mag- 
netic field, one can check that the expansions using a spin 
rotation-invariant formalism or a symmetry-broken one 
provide equivalent results for the physical observables (at 



5 = 0), such as the total T-matrix [7f 



IV. RESULTS 



In this section we investigate the analytic expressions 
obtained previously, both for antiferromagnetic and fer- 
romagnetic exchange, and compare them to numerical 
renormalization group (NRG) calculations'^ The 
simulations will be performed with the full density matrix 
using the DM-NRG algorithm^^''^'', which proved more 
accurate to compute observables in a finite (and possi- 
bly large) magnetic field compared to the standard NRG 
procedure. 



A. Anti-ferromagnetic Kondo model 

The weak coupling RPA result Eq. (38) for the T- 
matrix obtained at zero magnetic field can be rewritten 
in the limit D ^ \Ld\ ^ Tk in terms of the Kondo tem- 
perature Tk'- 



-upoT^ ^(w) 



37r^ 



1 



16 



loe 



(52) 



This recovers the standard poor man's scaling result, 
with some inclusion of high logarithmic corrections, 
which cut-off the divergence at uj — Tk (nevertheless, 
the RPA breaks down in this regime as well). This dif- 
ference between the two results stems from the fact that 
the imaginary part of the polarization diagram is not ne- 
glected for the renormalization of Kondo exchange in the 
RPA. A benchmarking of this expression to the DM-NRG 
calculation was already provided in Fig. 5. 

We now concentrate on the case of a finite magnetic 
field, where the presence of two scales, u and -B, prevent 
a simple poor man's scaling, and requires more sophis- 
ticated analytical renormalization group techniques^^'^^. 
Our result Eq. (50) obtained with the Majorana can be 
re-expressed in the regime D ^ u ^ B or —D <^ oj ^ B 
as: 



leiog^d^l) 



16V(|7^|) 



for D:^ uj:^ B 



for - £» < a; < B 



(53) 



resulting in different asymptotes in the case of large posi- 
tive and negative frequencies, in agreement with the fRG 
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FIG. 8. (color online) Spin-resolved T-matrix of the Kondo 
model with antiferromagnetic exchange in a finite magnetic 
field B = 5000Tk as computed in the Majorana diagrammat- 
ics at RPA level from Eq. (50). The high energy asymptotes 
given in Eq. (53) are also shown. A remanence of the Kondo 
peak occurs for frequencies near uj — B, see details in Fig. 9. 



results by Rosch et al}^'^^. A global view on the the T- 
matrix is given in Fig. 8, which shows a resonance that is 
now centered at uj = B. Note that perturbation theory 
is adequate because the scattering process remains far 
away from the unitary limit, —npoT ^{^) ^ 1- 

Let us now consider the analytical results in more de- 
tail, by comparing them to DM-NRG simulations. For 
all the NRG calculations, we take the half bandwidth 
D — 1 and the Kondo coupling J/D — 0.145, so that 
the Kondo temperature is Tk/D ~ 10^^. All results will 
be obtained using the U{1) symmetry and keeping the 
first 320 states, using a logarithmic discretization of the 
band with NRG parameter A = 2. We consider also the 
standard broadening parameter b in the range 0.5 to 0.9, 
in order to check possible problems of convergence in the 
broadening method (which can arise especially near the 
resonance at lu — B). Fig. 9 shows the comparison of the 
NRG calculation and the analytical expression Eq. (50) 
for the 7~-matrix. Clearly the agreement is overall excel- 
lent, but let us discuss rather the nature of the visible dis- 
crepancies. The most evident disagreement occurs near 
the resonance at oj — B, where the NRG shows a very 
broad shoulder while the Majorana result at RPA level 
presents a threshold followed by a spurious logarithmic 
divergence. The study of several broadening parameters 
in Fig. 9 shows that the broadening of the NRG data con- 
verges well except near the resonance, so that the actual 
result is likely to display a threshold-like feature at uj = B 
similar to the Majorana result. This broadening issue 
of the NRG is a notoriously difficult problem, to which 
more refined broadening methods^"'^'^^ do not bring much 
improvement^^'"*"^. In contrast to this dynamically gen- 
erated resonance of the T-matrix, atomic like resonances 
(which occur in the abssence of dissipative coupling to 
the electronic environment) can be tackled by improved 
broadening method'*'^^''*^, as we will see in Sec. IV C for 
the transverse spin susceptibility. Despite this flaw of the 
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FIG. 9. (color online) Close up on the finite magnetic field res- 
onance in the Kondo model with antiferromagnetic exchange. 
This is a comparison between Majorana diagrammatics at 
RPA level (dashed line) and DM-NRG simulations with de- 
creasing broadening parameter b — 0.85, 0.75, 0.65, 0.55 (solid 
lines, top to bottom) for a ratio B/Tk = 5000. The NRG data 
are clearly overbroadened for frequencies near the threshold 
u) = B. One note also a slight overshoot of the Majorana re- 
sult at low frequency, also visible in Fig. 8, an artifact which 
is discussed in the text. In all other frequency ranges, the 
agreement between analytics and numerics is quantitative. 



NRG to predict a threshold aX uj — B, the NRG simula- 
tion predict correctly a finite amplitude peak above the 
threshold, while the Majorana diagrammatics leads to a 
spurious logarithmic divergence at w — > _B+, which can 
be clearly traced back to the form of the /-fermion bubble 
in Eq. (46). It has been argued in the context of fRG^^'^^ 
that this spurious singularity should be smoothened by 
taking into account spin relaxation effects, an issue that 
we examine in Sec. IV C. 

A second drawback of the present RPA resummation 
scheme can be observed in the low frequency range in 
Figs. 8 and 9. While the plateau of the T-matrix in the 
range Tr- ^ w ^ -B is quantitatively reproduced by the 
Majorana diagrammatics, a spurious logarithmic diver- 
gence occurs at uj <^ T^:, which leads to progressive de- 
viations from the numerical results at lowering frequency. 
The origin of this behavior is easily understood from the 
Majorana fermion bubble Eq. (45), which displays an 
unbounded logarithm in the zero frequency limit. We 
have checked however that including spin relaxation ef- 
fects will not cure this divergence, which has rather to 
do with deficiencies of the RPA resummation. Indeed, 
while the RPA diagrams reproduce exactly the one-loop 
poor man's scaling approach at zero magnetic field, the 
RPA scheme at finite magnetic field is actually not fully 
equivalent to the fRG approach in the limit uj <C Tk 
(the two approaches give however equivalent results oth- 
erwise). This difference can be understood from the fact 
that a one-loop renormalization of the Kondo exchange 
takes into accounts non-RPA diagram, such as seen in 
Fig. IV A. One can check that these non-RPA diagram 
remain logarithmic in the range w ^ Tk for finite mag- 
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FIG. 10. Upper panel: standard one-loop renormalization of 
the Kondo exchange. Resummation of the top right diagram 
(particle-hole channel) can be generalized in the Majorana 
diagrammatics in terms of the non-RPA extra contribution 
shown in the lower panel. The neglect of such terms which 
remain logarithmic in the range uj <^ Tk for finite magnetic 
field is the likely origin for the discrepancies of the RPA at 
low frequency (and finite B) with respect to one-loop fRG 
calculations. 
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FIG. 11. (color online) Kondo resonance in a finite magnetic 
field for antiferromagnetic exchange, with the same param- 
eters as Fig. 9. This plot shows the corrected RPA with 
the approximate particle-hole channel (RPA+p-h) diagram 
of Eq. (54), with comparison to NRG calculations and the 
uncorrected RPA. 



netic field, and these will likely cancel out the spurious 
low frequency divergence of the RPA scheme (these dia- 
grams are however logarithmically subdominant at high 
energy, allowing the good agreement between RPA and 
fRG for the other frequency regimes). We have unfortu- 
nately not been able to resum exactly to all orders the 
non-RPA diagram shown in the lower panel of Fig. IV A. 
However, it is possible to put the present argumentation 
on a more precise ground, by noting that the leading con- 
tribution of the shown diagram in Fig. IV A is given by 
the combination [no,,(w)]^no/(w — 0). Iterating these 
terms in a geometric series to all orders amount to can- 
cel the remaining spurious logarithmic divergence at low 
energy, as given by this new form of the T-matrix: 



T';"ito) 



J2 

— 8~ ^^S^*-^-* ^ 



(54) 



Im 



n, 



2n, 



0/ 



1 



n, 



0/ 



^ [nojjiio/ 



n;„no/(o)]. 



Figure 11 shows indeed the elimination of the spurious 
low-w resonance, and quantitative agreement with the 
NRG for all frequencies (except right at the threshold 
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FIG. 12. (color online) T-matrix of the Kondo model with 
ferromagnetic exchange J/D — —0.138 and finite magnetic 
field B/D = 0.001, comparing the DM-NRG calculations for 
broadening parameters b = 0.85,0.75,0.65,0.55 (top to bot- 
tom) to the RPA result with particle-hole corrections Eq. (54). 

bj = B). This successful benchmarking of the Majorana 
diagrammatics allows us to consider now with confidence 
the T-matrix in the case of ferromagnetic exchange. 



B. Ferromagnetic Kondo model 

We now turn to the central question of the present 
work, namely the magnetoconductance of the spin S = 
1/2 Kondo model with ferromagnetic exchange. This de- 
scription also apphes to the strong coupling hmit of the 
spin S = 1 underscreened Kondo model with antiferro- 
magnetic exchange^^'**^, recently observed in molecular 
quantum dots'^^^^'^'*. For both models, the ground state 
remains a spin doublet with logarithmic decoupling from 
the bath of conduction electrons, and so one expects an 
extreme sensitivity to an applied magnetic field. The 
ferromagnetic exchange modifies dramatically the flow 
equation of the Kondo coupling, which remains small at 
all energies, and even renormalizes to zero in the low fre- 
quency limit 1*5,44^ -y^g compare again the RPA results 
(with particle-hole corrections) with DM-NRG simula- 
tions, see Fig. 12. We find excellent agreement between 
the numerics and Majorana perturbation theory, as ex- 
pected from the perturbative nature of the ferromagnetic 
Kondo model. However, near the resonance at w = S, 
one observes again that the DM-NRG results are over- 
broadened, similar to the antiferromagnetic case, so that 
the step predicted by the theory is difficult to repro- 
duce. Perturbation theory displays also an artifact, as 
it predicts that the T-matrix vanishes logarithmically at 
Lo — a, behavior that will be corrected by including 
spin relaxation effects, see Sec. IV C. 

We finally address the issue of the Zeeman splitting in 
the spin-symmetric T-matrix, see Fig. 13 which compares 
experimental results for a = 1 underscreened molecular 
quantum dot^^ and our analytical result. The analytical 
formula shows indeed that the Zeeman splitting persists 
for arbitrary small values of the applied magnetic field. 
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FIG. 13. (color online) Upper panel: experimental data^'^ for 
the finite bias conductance of an underscreened 5* = 1 Kondo 
dot as a function of magnetic field and applied bias. A split- 
ting is still observed at small magnetic field, below the typical 
voltage scale characterizing the width of the B = differ- 
ential conductance. Lower panel: spin-symmetric T-matrix 
of the 5* = 1/2 Kondo model with ferromagnetic exchange 
{J/D — —0.138) as a function of magnetic field and frequency. 
Similarly to the experiment, the Zeeman effect persists even 
for vanishingly small magnetic field. 



giving strength to the interpretation made in Ref. 13. 



C. Spin relaxation time 

We finally investigate the role of spin relaxation in 
transport, motivated by the spurious logarithmic singu- 
larities obtained at RPA level in the T-matrix near the 
threshold uj = B. Following previous works^*"^^, one 
can expect that spin lifetime effects occuring at higher 
order in perturbation will smear out the artificial diver- 
gence of the T-matrix. In the Majorana diagrammatics, 
spin relaxation in straightforwardly incorporated within 
Majorana self-energies, because the magnetic susceptibil- 
ities identically reduce to Majorana propagators, see e.g. 
Eq. (20). 

Due to the presence of the magnetic field along the 
z-axis, we need to introduce longitudinal and trans- 
verse susceptibilities, Xz{t) = — {Tt[S^{t)S^{Q)]) and 
X_l(t) = -{Tr[S+{T)S-{Q)]). The spin operators are 
then replaced by = $77, S+ = V^^P and S~ = 
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\/2$/, so that imaginary parts of these susceptibilities 
are given by single fermion propagators: 



X?(a;) = -itanh(^)G;«H (55) 

X?(w) =-tanh(^)G;«(a;). (56) 

Spin relaxation is now immediately incorporated from 
Dyson's equation: 



1 

OJ — Sr)(w) + iO+ 

1 



(57) 
(58) 



and the needed imaginary part of the sclf-cncrgics arc 
easily computed up to second order in J (here at zero 
temperature): 



I];«(a;) = 



J2 



-Po 



U'^fiuji +oj) - liofi-uji - w) [Sign(wi + ui) - Sign(wi)] dwi 



t"R 



I rD 

-Po / U'ofiu 

J-D L 



wi + w) + Uq:^{u}i + w) [Sign(a;i + w) - Sign(wi)] dcoi, 



(59) 
(60) 



from which the real part of the self-energies can be ob- 
tained by applying the Kramers-Kronig relation. The 
self-energy contains an extra w-independent contri- 
bution given by 

SEr, = dTGoc{T)Goc{-T) = ^ po log(4), (61) 

that originates from the last term in Eq. (41), and which 
is crucial to preserve the correct particle-hole symmetry. 
The susceptibilities obtained from the Majorana propa- 
gators (55)-(56) are shown in Figure 14. 

Let us begin with the discussion of the transverse s\is- 
ceptibility, which is expected to display a narrow res- 
onance at the Zeeman energy. As noted previously^'', 
the NRG data for the transverse spin susceptibility is 
much more overbroadened than the T-matrix (which only 
shows a shoulder and not a peak), as can be gathered 
from the sharper feature found in the Majorana calcu- 
lation. Some improvement of the NRG data processing 
can be obtained by using adaptative broadening methods 
("b-trick")^^, where the broadening parameter b is taken 
to be freqiiency-dependent. This allows to obtain a much 
narrower peak, yet still somewhat too broad. Larger scale 
NRG simulations using so-called 2;- averaging^ ^ should al- 
low to obtain better agreement. We now turn to the lon- 
gitudinal susceptibility, which similarly to the T-matrix 
is expected^^ to give a shoulder atu} = B. The lowest or- 
der Majorana calculation shows good agreement with the 
NRG at a; > S, but fails to reproduce the low- frequency 
tails, because the self-energy Eq. (60) has a gap. Previous 
work^"^ showed however that enforcing self-consistency 
in the perturbative scheme allows to cure this defect of 
the bare perturbation theory. Note also that the use of 
self-consistency in the self-energies allows to include the 



rcnormalization of the exchange coupling within the re- 
laxation rate, a small quantitative effect that we did not 
take into account here. 

Ultimately, it is the finite amplitude of the resonance 
peak in the spin susceptibilities at uj — b which is ex- 
pected to regularize the spurious divergence in the T- 
matrix^®"^^. This is simply achieved by replacing the free 
propagators Eq. (43) by the dressed ones Eq. (57)-(58) 
into the previous T-niatrix calculation. Figure 15 shows 
the inclusion of spin relaxation effects at the RPA level 
for the T-matrix for the ferromagnetic Kondo model. 
The previously observed artifact, namely a logarithmic 
cancellation of T ■'^ at w = S, is now cured by the spin 
lifetime. Note still that the shoulder at the Zeeman en- 
ergy persists, an important check to argue about the oc- 
currence of a split Kondo resonance at low magnetic field 
in the underscreened situation^^. 



V. CONCLUSION 

The goal of this paper was to imdcrstand quantum 
transport under an applied magnetic field in presence 
of a ferromagnetic Kondo interaction. Although mag- 
netic impurities (or spin-active quantum dots) are most 
often coupled antiferromagnetically to conducting elec- 
trons (see however Rcf. 34), the ferromagnetic case is 
relevant for the situation of imdcrscreened large mag- 
netic moments^^'^*. Indeed, partially screened moments 
decouple weakly from the electrons at low energy, due to 
an effective ferromagnetic Kondo interaction^^. The ex- 
istence of a degenerate magnetic ground state leads then 
to a singular response to magnetic field, also visible in 
transport, which we investigated here in some detail, us- 
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FIG. 14. (color online) Upper panel: transverse spin suscep- 
tibility for the antiferromagnetic Kondo model, comparing 
the lowest order Majorana calculation to NRG simulations 
obtained with constant broadening"^"^ b = 0.55 and with the 
adaptative broadening methods ("b-trick")*^, using the same 
parameters as in Fig. 9. The resonance at oj = B is clearly 
overbroadened, and the broadening trick allows to improve 
on the resolution. Lower panel: longitudinal spin susceptibil- 
ity within the Majorana diagrammatics and the NRG (with 
constant broadening b = 0.55 only). The spurious gap ob- 
tained by the lowest order diagrammatics result should be 
filled by enforcing self-consistency in the diagrammatic calcu- 
lation scheme'^^. 



posed various renormalization scheme, either based on 
Callan-Symanzik equations'*^, functional RG^^'^^, flow 
equations^*^ , and Liouvillian-based RG^*'^^. However 
we have to recognize that the Majorana diagrammat- 
ics is likely harder to implement beyond one-loop, be- 
cause disentangling the sub-leading logarithmic contri- 
butions really requires an RG method (the Majorana 
scheme performs yet very well beyond one loop near a 
perturbatively-accessible quantum critical point^''). An- 
other issue addressed here and in previous works**'"'^^ is 
the energy resolution of the NRG, which is clearly in- 
sufficient. We believe that a careful comparison of large 
scale NRG simulations*^'"'^'*^ and possibly two-loop RG 
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FIG. 15. (color online) T-matrix of the ferromagnetic Kondo 
model near the Zeeman energy B = 0.005-D, comparing the 
RPA calculation with and without self-energy corrections. 
The inclusion of a finite spin relaxation rate allows to cut 
the singularity (a spurious logarithmic vanishing at oj = B) 
of the bare Majorana calculation. 



analytics^"'^'^^, already in the equilibrium situation, still 
need to be performed. 
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